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Abstract 

We investigate the one-gluon-exchange (acts) corrections to the po- 
larized real photon structure function gi{x,Q'^) in the massive parton 
model. We employ a technique based on the Cutkosky rules and the 
reduction of Feynman integrals to master integrals. The NLO contri- 
bution is noticeable at large x and does not vanish at the threshold of 
the massive quark pair production due to the Coulomb singularity. It is 
found that the first moment sum rule of gj is satisfied up to the NLO. 
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The experiments at the Large Hadron Colhder (LHC) have started and it is much 
anticipated that signals for the Higgs boson and also for the new physics beyond the 
Standard Model (SM) will be discovered [1]. Once these signals are observed, more 
precise measurements will need to be performed at the future e+e^ collider, so-called 
the International Linear Collider (ILC) [2]. In such detailed knowledge of 

the SM at high energies, especially based on QCD, is still important. 

It is well known that, in high energy e~^e~ collision experiments, the cross section 
of the two-photon processes e~^e~ — )■ e'^e" +hadrons dominates over other processes 
such as the annihilation process e~^e~ — t- 7* — hadrons. The two-photon processes 
at high energies provide a good testing ground for studying the predictions of QCD. 
In particular, the two-photon process in which one of the virtual photon is very far off 
shell (large = —q^), while the other is close to the mass shell (small = —p^), 
can be viewed as a deep-inelastic electron-photon scattering where the target is a 
photon rather than a nucleon. In this deep-inelastic scattering off a photon target, 
we can study the photon structure functions, which are the analogs of the nucleon 
structure functions. When polarized beams are used in e^e~ collision experiments, 
we can get information on the spin structure of the photon. 

For a real photon (P^ = 0) target, there exists only one spin-dependent structure 
function gl^x, Q^), where x = Q"^/ {2p-q). The photon structure functions are defined 
in the lowest order of the QED coupling constant a = e'^/An and they are of order 
a. The QCD analysis of gj was performed in the leading order (LO) (the order a) 
[3], and in the next-to-leading order (NLO) (the order aa^) where = g^/4:TT 
is the QCD coupling constant. In these analyses all the active quarks are treated 
as massless. At high energies the heavy charm and bottom quarks also contribute 
to the photon structure functions. The NLO QCD corrections due to heavy quarks 
have been calculated for the unpolarized photon structure functions F2 {x, Q^) and 
F2{x, Q^) [5j. The heavy quark mass effects on gj were analysed at NLO in QCD in 
Ref.[6] by using the LO result of the massive parton model (PM). But the complete 
heavy quark mass effects have not yet been computed for gj at NLO. 

In this paper we investigate the real photon structure function gJ in the massive 
PM at NLO in QCD. In order to compute gJ at NLO, we employ a technique based 
on the Cutkosky rules [7] and the reduction of Feynman integrals to master integrals. 
The master integrals which appear in this analysis also show up in computing other 
photon structure functions such as F2{x,Q^) and F2{x,Q^) at NLO. We express 
the phase space integrals of these master integrals in analytical form as much as 
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possible so that they may serve as useful tools for the analyses of the future ILC 
physics. 

The polarized real photon structure function gl satisfies a remarkable sum rule 

i El [101 [m [12] 

' g2{x,Q'^)dx = . (1) 

In particular, applying the Drell-Hearn-Gerasimov sum rule [13] to the case of a vir- 
tual photon target and using the fact that the photon has zero anomalous magnetic 
moment, the authors of Ref. [T2| argue that the sum rule ([1]) holds to all orders in 
perturbation theory in both QED and QCD. We examine whether the NLO result of 
gj in the massive PM satisfies this sum rule. We find numerically that the sum rule 
([1]) is indeed satisfied at this order. But we point out that the sum rule may not be 
well-defined when gj is analysed to higher orders in perturbation theory, since the 
calculated result may diverge at the threshold of the massive quark pair production 
due to the Coulomb singularity. 

We calculate the cross sections for the two photon annihilation to the heavy 
quark gnQH pairs 

i*{q) + i{p) qH + QH , (2) 

with one-loop gluon corrections and to the gluon bremsstrahlung processes 

7*(g)+7(p) — >qH + qH + 9- (3) 



We employ the technique developed by Anastasiou and Melnikov [M], which is based 
on the Cutkosky rules and the reduction of Feynman integrals to master integrals. 
First, following the Cutkosky rules p], the delta-functions which appear in the phase 
space integrals are replaced with differences of two propagators 

2m6{r' - m') ^ - , (4) 

where m is the heavy quark mass. Then the cross sections for the virtual corrections 
to the processes ([2]) and for the bremsstrahlung processes ([3]) are described by the 
two-loop diagrams shown in Fig{T] and FigJJl respectively, where a cut propagator 
should be understood as the r.h.s. of Eq.(|l]). 

We regularize the amplitudes by dimensional regularization D = 4 — 2e. Then 
we apply the following D-dimensional projection operator 
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Figure 1: Two- loop diagrams with virtual corrections. Graphs with virtual correc- 
tions to the right of the cut lines and graphs with (g, /i) and {p, p) interchanged are 
added. Graphs with the external quark self-energies are not shown in the Figure, 
but should be included in the calculation. 
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Figure 2: Two-loop diagrams with a real gluon emission. Similar graphs correspond- 
ing to (e) and (f) are included. Also graphs with (g, /i) and p) interchanged are 
added. 
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with 

to these diagrams to extract the contributions to gj. They are expressed in terms 
of a large number of two-loop scalar integrals of the form 

A{Vi) = A{Vk, ^kq, ^kp, ^kpq, ^l, J^lq, J^lp, ^Ipq, l^kl) 

d^k d^l 1 

(27r)^ (27r)« - w?^' [{k - qf - m?]"^-^ [{k - pf - m^Y'^^ [{k - p - qf - 

1 



X 



[P - m^Y' [(/ - g)2 - m?]"'" [{I - pf - \{l - P - qf - m^^'^ \{k - OT" ' 

(7) 

and the coefficients of these integrals are written as functions of a;, Q^, rn? and D. 
Note that Xjik — is a gluon propagator. Actually Aiv^ has seven propagators at 
most and thus at least two z/j's are zero. 

We arrange the integration variables k and / so that the cut propagators are 
— m^] and — p — g)^ — m^] for the diagrams in FigJH Among many v4(z/j)s 
there appear those with one or both of the cut propagators eliminated. Those 
integrals do not contribute to g{. Thus we only pick up Aiv^'s, which are in the 
form z/feq, i/fcp, 1, z/;, i/jq, i/jp, i/^p,, i/fcj) and discard others. A similar procedure is 
applied to the diagrams in FigJ2J We choose — vp?\^ ~ P ~ qY ~ ^'^] 

and l/{k — for the cut propagators and, therefore, search A(z/j)s in the form 
A{i>k, ^kq, ^kp, 1, 1, ^iq, ^ip, ^ipq, 1) aud dlscard others. 

The number of the relevant A(z/j)s is still large. Then, following the reduction 
procedure [l5] which is based on the method of integration by parts [T6] and the use 
of the Lorentz invariance of scalar integrals [17], these v4(z/j)s can be expressed in 
terms of fewer number of master integrals. Today several public codes fTEj [T^ I2U] 
are available. We make use of FIRE and express the relevant A{i>i)s as a linear 
combination of the master integrals which are denoted as 

M{Ui) = M{Uk, Vkq, l^kp, l^kpq, l^l, l^lq, l^lp, l^lpq, m) , (8) 

in the same way as the notation of A(z/j)s in Eq.(I7]). Again the master integrals in the 
form of M(l, u^q, ^kp, 1, ^iq, ^ip, ^ipqj ^ki) are only relevant for the virtual correction 
diagrams in FigH] and those in the form of M(i/fc, i/^g, i/^p, 1, 1, i/^g, z//p, i/^pg, 1) are 
relevant for the real gluon emission diagrams in Figj2l 
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Finally we perform the phase space integrations for these cut master integrals. 
For the two-cut and three-cut master integrals, we evaluate 



DisC^^) M(l, Z/fcq, I/fcp, 1, Vl, Uiq, Uip, Uipg, Vkl) 

d^l 1 



(9) 

and 

DisC^^) M(z/fc, Z/fcq, Z/fcp, 1, 1, Z//q, Z/;p, Z//p,, 1) 



LTX] 

1 



X 



- m^'' \{k - qf - m^'"' \{k - - w?]"''^ 

1 

^ [(/ - g)2 - m?]"'" [{I - pf - m^'^ \{l - P - qY - m^l"'''^ ' 



respectively. Note that at least two z/j's are zero in both ([9]) and (UOl) . 

There appear 61 master integrals in total in this analysis of g[[x^ Q"^). However, 
the choice of a set of master integrals is not unique. We are at liberty to replace 
a master integral with one of the other scalar integrals. We choose a set of master 
integrals such that each corresponding coefficient function is finite in the limit D — )■ 
4 [21]. With this choice of the set, the phase space integrations for master integrals 
need only be evaluated up to the finite terms in the series expansion in e. 

When the virtual correction diagrams in FigJT]are concerned, the ultraviolet (UV) 
singularities appear in the graphs (b), (c) and (d), while the infrared (IR) singularity 
emerges from the graph (a). Both the UV and IR singularities are regularized by 
dimensional regularization. The UV singularities are removed by renormalization. 
We adopt the on-shell scheme both for the wave function renormalization of the 
external quark and the mass renormalization. For the latter, we replace the bare 
mass in the Born cross section by the renormalized mass m. 



Air \ / I e 



(11) 



where Cp = | is the Casimir factor, S'' = (47r)^e ^"'^ with Euler constant ■je and /i 
is the arbitrary reference scale of dimensional regularization. The renormalization of 
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Figure 3: Charm quark effects on the polarized real photon structure function 
c/7(a;,Q2) in the PM in units of (3ae^/7r) for = 20 GeV^ and = 1.3 GeV 
with as = 0.22. We plot the LO result (dotted line), the NLO contribution (dashed 
line) and the sum of LO and NLO contributions (solid line). 



the QCD gauge coupling constant is not necessary at this order. The IR singularities 
appear also in the real gluon emission graphs (a),(b), (c) and (d) of Figj2l However, 
the IR singularities cancel when the both contributions from the virtual correction 
graphs and the real gluon emission graphs are added. Actually the IR singularities 
reside in the two-cut master integrals in the form M(l, u^q, ^kp, 1, 1, ^iq, ^ip, 1, 1) and 
the three-cut master integrals Miy^, u^g, I'kp, 1, 1, ^iq, ^ip, ^ipq, 1) with z/^ + uipg = 2. 
The details of the calculation will be reported elsewhere j,22j . 

In Figs I3] and m we plot the polarized real photon structure function gl^x^Q"^) 
predicted by the massive PM up to the NLO for the case of = 20 GeV^ and 
tts = 0.22. We choose c and 6 as a heavy quark, for Figs. [3] and HI respectively. We 
take mc = 1.3 GeV, rrij, = 4.5 GeV, Cc = | and Cb = — |. Here we show three curves: 
the LO result, the sum of LO and NLO corrections and the NLO corrections alone. 
The allowed x region is < x < a;niax with 

_ 1 

The LO result is expressed by 

^77(^,Q')|LO = ^e:|(ln[^)(2x-l)+/^(-4x + 3)| , (13) 
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Figure 4: Bottom quark effects on the polarized real photon structure function 
g1{x, Q2) in the PM in units of (Sae^/vr) for = 20 GeV^ and mb = 4.5 GeV with 
as = 0.22. We plot the LO result (dotted line), the NLO contribution (dashed line) 
and the sum of LO and NLO contributions (solid line). 



where 

For X — )■ Xmax, goes to zero and thus (?! (x, (5^)lo vanishes at Xmax- 

We observe in the Figures that there exist NLO corrections both at large and 
small x, positive at large x but negative in small x region, a behavior similar to the 
LO result. Especially, the radiative corrections are large near the threshold (near 
Xjnax) and the NLO curve does not vanish at a;max- This is due to the well-known 
Coulomb singularity, which appears when the Coulomb gluon is exchanged between 
the quark and anti-quark pair near threshold. The diagram FigHJ^a) is responsible 
for this threshold behavior. The virtual correction to the left of the cut line in 
Fig{T]^a) gives rise to a factor while a factor /3 comes out from the phase space 
integration. They are combined and yield a finite but non-zero result at Xmax- 

We consider the sum rule ([T]) for a real photon target. Substituting the LO result 
i/i (a;, Q^)|lo given by ([13]), we see the sum rule holds [H El HOl [IH [I2] . Figsj3]and|l 
show that the sum rule also seems to be satisfied by the NLO contribution to g{ in 
both cases. Expressing the NLO contribution as (71 (x, (5^)|nlo! "we find numerically 

/ gJix,Q^)\^i^odx = . (15) 
Jo 
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But due to the limitation of accuracy of our numerical integration, we observed 
^ ^ k"-!>l(-.Q')Wod. ^ (_2.2, - 2,0) X 10- . 



NLO 



dx 



for charm and bottom cases, respectively. So we conclude that the sum rule is 
satisfied in the massive PM up to the NLO. 

However, if we go on further and analyse gi{x,Q'^) to higher orders in pertur- 
bation theory, we expect that the result will diverge at Xmax due to the Coulomb 
singularity. A detail analysis on the structure of the Coulomb singularity tells that 
S'l Innlo ~ /3 X [as/ fi)^ [231 121] whose integral for the first moment is ill-defined 
due to end-point singularity at x = Xmax- The sum rule is not well-defined in the 
perturbation theory starting at NNLO. To obtain an appropriate threshold behavior 
for photon structure functions, we may resort to the method of resummation of the 
Coulomb singularities. A noticeable difference in the resummation is emergence of 
bound-state poles of qu^H above Xmax- Then the left-hand side of the sum rule Eq. 
([1]) should include also the bound-state contributions. We will not pursue this issue 
further here but render it to our future publications. 

In summary we have calculated the NLO corrections to the polarized photon 
structure function gl in the massive PM. We have found that the NLO contribution 
is noticeable at large x and does not vanish at Xmax due to the Coulomb singularity. 
We have also found numerically that the sum rule ([1]) is satisfied up to the NLO in 
the massive PM. The details of our calculation will be reported elsewhere |22] . 
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